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ABSTRACT 


Monterey Bay, on the west coast of the United States, is unique in 
that it is a large, symmetric, semi-elliptical bay divided north and 
south by the deep Monterey Canyon. The effect of the canyon on seiching 
within the bay and on long wave oscillations within the bay was studied 
by anelleniae synchronized wave records at each end of the bay. Power 
Spectra and cross spectra calculated for five periods selected from six 
months continuous data indicate the Monterey Canyon has a profound effect 
on the bay's oscillating characteristics. The canyon appears to act as 
an impedence barrier dividing the bay into two independent oscillating 
basins each having recurring long-period waves which persist during 


significant long-wave activity. 
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I. INTRODUCTION 


Monterey Bay (Figure 1) lies on the west coast of the United States, 
located about sixty nautical miles south of San Francisco. The Bay is 
a large semi-elliptical bay which has some topographical features which 
are unique. A deep trough of the North Pacific basin, bounded on the 
north by the Mendocino Seascarp and on the south by the Murray Sea- 
scarp, approaches closer to the et at Monterey Bay than at any other 
point along the North American coastline. As a result, long waves propa- 
gated across the Pacific Ocean in an easterly direction may be somewhat 
contained. The continental shelf within the bay extends out to approxi- 
mately the 600 foot depth contour and is cut by the deep Monterey Canyon 
which has a volume of over 50 cubic miles. Monterey Canyon will be seen 
to exert significant influences on the oscillating characteristics of 
the Bay. 

There have been several studies concerned with the effects of long- 
period waves in Monterey Bay (a long-period wave for this study is 
defined as a wave having a period in excess of 1 minute), but these 
Studies have not been able to define the causes of the long wave 
activity and have left many areas of inquiry. Hudson [1947] proposed 
a surge model of Monterey Harbor for the U. S. Army Corps of Engineers. 
He used six months data (October 1946 - April 1947) from three elec- 
tronically synchronized sea-surface recorders. These results are 
tabulated in Table I. The presence of the long-period waves was 
observed. However, the cause of the surge mechanism was unknown. 

Wilson [1965] conducted a long-wave study of Monterey Bay for 


the U. S. Army Corps of Engineers to determine the feasibility of 
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construction of an engineering model of the surge phenomenon that 
occurred at various times within Monterey Bay. Statistical data were 
aaroe” from sensors located at Monterey Harbor, and Santa Cruz 
feecure 1). Three long-period wave recorders were also installed 
within Monterey Harbor and were in continuous operation from October 
1963 to April 1964. Two were arranged so that tides and sea-swell 
were filtered out, and a third recorded swell approaching the harbor. 
Wilson concluded that the surge within Monterey Harbor was not a result 
of incoming swell as no correlation was found to exist between long- 
period waves within Monterey Harbor and the incoming sea-swell. The 
Santa Cruz sensor was designed so that it functioned as a long-wave 
recorder. The recorder operated continuously from October 1963 to 
February 1964. Wilson's results for Monterey Harbor are summarized 
in Tables II and III. 

Wilson performed a residuation analyses of different records from 
sensors located around the bay in order to determine a local evaluation 
of the oscillations within the bay. Residuation analyses is accomplished 
by successively subtracting "apparent" periods from the wave record 
until a smooth trace remains. The subjectivity of which this method is 
accomplished can lead to ficticious results. These results are listed 
in Table IV. Table V is a synopsis of spectrum analyses for three days 
record of the Monterey sensors. Wilson concluded, based on Tables IV 
and V along with calculated modes of oscillation, Figures 24-29, that 
Monterey Canyon functions as an impedence between north and south 
portions of the bay, and therefore, free oscillations are to a large 


extent uncoupled. Wilson further states that the effect of any sharp 


discontinuity in submarine topography, equivalent to the edge of a 
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continental shelf, is to serve as a nodal position for any shelf 
oscillation to which the shelf is susceptible. 

(Robinson [1969] analyzed long-wave activity in Monterey Bay on 23 
Ronuaty and 20 April 1969. Tide recorders were located at Monterey 
Wharf #2 and Santa Cruz Municipal Wharf, however, the records were not 
synchronized. Thus, cross spectra and phase information was not 
obtainable. Robinson concluded that Monterey Canyon affected the 
seiching motions within the bay, however, it did not appear to divide 
the bay into two independent oscillating basins, as concluded by Wilson 
[1965]. Similar periods and amplitudes were found at Monterey and 
Santa Cruz concurrently, but correlation between the two locations 
could not be computed due to non-synchronization of records. 

In the present study simultaneous tidal records from tide gages 
located at Monterey Wharf #2 and Santa Cruz Mucicipal Wharf were 
analyzed during periods of significant long-wave activity in Monterey 
Bay. Continuous analog records were maintained from 3 November 1969 
to 30 April 1970 and five days (6 November, 3 December, 6 December, 
14 January, 20 January) were selected for analysis and comparison. 

The objectives were to perform an overall long-wave study of 
Monterey Bay, and to examine various modes of free oscillation of the 
bay to determine if in fact seiching does exist in the bay. This was 
accomplished through examination of the individual power spectra at 
Monterey and Santa Cruz, computation of the cross spectra and deter- 


mination of phase differences of long-period waves present. 
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TABLE I 


Long-Period Waves in Monterey Harbor (after Hudson, 1949) 


PERIOD AVERAGE HEIGHT PERCENTAGE OF 
(min) (ft) TIME PRESENT 


0.4 
O25 


not given 





TABLE II 


Marine Advisor's Data for Monterey Harbor (after Wilson, 1965) 





SENSOR PERIOD AVERAGE HEIGHT PERCENT OF 
(min) (ft) TIME PRESENT 
1 iL Foil 0.1-2.5 0-50 
2 1-14 0.1-3.0 0-55 
TABLE III 


Marine Advisor's Data for Santa Cruz Harbor (after Wilson, 1965) 


PERIOD AVERAGE HEIGHT PERCENT OF 
(min) (ft) TIME PRESENT 
e203 0-18 
2.3-4 0.2-2.0 20-70 
10-14 0-25 
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II. INSTRUMENTATION 


A. MONTEREY 

The data for Monterey was obtained using a standard Coast and 
Geodetic Survey automatic tide gage. (Manual of Tide Observations, 
1965), located on Monterey Municipal Wharf #2. Monterey Harbor. The 
tide gage is maintained daily by NPS personnel. 

The Monterey tide gage senses changes in water level by means of 
a float/pulley arrangement, which is entirely mechanical. The recording 
drum is advanced by a clock mechanism. The drum speed is designed for 
1 in/hr, however, the hourly feed was measured to be 1.03 in/hr. The 
instrument is time checked daily with time marks being accurate within 
+ 5 seconds. The marigram is recorded in rectilinear coordinates on 
plain white paper. A stilling well, which is a 12 inch diameter steel 
pipe with a 1 inch orfice in the bottom, serves as low pass filter for 


the majority of wind waves (periods of one minute and below). 


B. SANTA CRUZ 

The Santa Cruz data were obtained with a Bristol Model 28 gas- 
purging pressure (bubbler), portable tide gage, (Figure 2), located 
on Santa Cruz Municipal Wharf. This instrument senses changes in 
water level by means of a nitrogen-filled tube which is connected to 
a bellows system (Manual of Tide Observation, 1965). 

A bubbler orfice chamber was connected to the end of the sensing 
tube to reduce sensitivity to short period wave action. There is also 
a bellows inlet needle valve which could be throttled to filter add- 


itional wave action from the record. 


18 








\ 


ln 


ee ew eee ae 
e 


Ay 


4 ah bh 


g 


a oie AL ee 


i — 
ny 


S) 





-OOOMONDALKRWNM> 









Bs y’ = 
: Ba Fee ta! by Fy, 
ie oot ogdl if [ | QO é “TF? 





NITROGEN BOTTLE 

PRESSURE REDUCING VALVE 

PRESSURE DIFFERENTIAL REGULATOR 
ADJUSTABLE NEEDLE VALVE 
TRANSPARENT BUBBLER CHAMBER 
PULSATION DAMPER 

STRIP CHART RECORDER WITH TRANSDUCER 
TELEMETERING PRESSURE TRANSMITTER 
TUBING 

BUBBLER ORIFICE CHAMBER 

STILLING WELL 


FIGURE 2 THE BUBBLER TIDE GAGE. 
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The pressure reduction valve is regulated to allow a constant 
pressure greater than that due to the maximum head anticipated over 
the orfice. The pressure differential regulator together with the 
shunt line provide for a constant pressure difference of about 3 
pounds per square inch across the flow regulator. Thus, the rate of 
flow can be adjusted to a desired constant value and is not dependent 
on the tide stage. Flow is through a silicone oil-filled bowl so that 
it can be monitored visually. A rate of 30-60 bubbles per minute was 
used for this study. 

Advantages of this scheme are that the sensitive elements of the 
gage need not be designed to withstand the underwater environment; and 
fluctuations in barometric pressure are not directly reflected on the 
record. 

A significant portion of the system capacity is at the submerged 
orfice, thus providing some high frequency filtering and protecting 
the small diameter supply line and orfice from marine fouling. 

A Bristol Company clock recorder provides an analog record. The 
record is on a 6 inch (15.2 cm) wide strip chart and is recorded in 
curvilinear coordinates. Paper advance is adjustable and 6 in/hr was 
utilized. The clock drive has an eight-day spring, and the record is 
sufficient for about seven days. Manufacturer's claims are that 
hysteresis and/or non-linearity are limited to 1% of the full-scale 
value. For sea water use, the instrument is calibrated to correspond 
to a specific gravity of 1.025. 

In addition to the above high-frequency filtering mechanisms, a 
stilling well was designed and installed to attain the desired 


filtering of high frequency and waves [Robinson, 1966]. The well was 
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constructed of a 20 foot section of polyvinyl chloride (pvc) 6 inch 
inside diameter pipe. It was capped on 16, 1/4 inch inside diameter 
holes drilled in the side. Copper sleeves were inserted to eliminate 
Puliae. The 16 holes provided the capability of increasing the orfice 
from a 1/4 inch to a 1 inch diameter opening. 

The response characteristics of the well were determined theoretic- 
ally for two orfice sizes and three different wave frequencies using 


the equation for the rate of water rise in a well [Robinson, 1969]. 


— : a 
dh ,/dt 0.6 a/A /2g (ho h.) 


where, 
a = orfice area 
A = well area 
0.6 = empirical orfice flow coefficient 
g = acceleration due to gravity 
h, = water height outside the well, i.e., the forcing function 


h, 
i 


water height insdie the well. 


The forcing function, h,, was chosen to be a simple sine function 


0? 
of unit amplitude, and of frequency equal to the wave frequency of 
interest. The initial conditions were h, = 0 att =0. The results 
are summarized in Table VI. 

The response characteristics for the Monterey stilling well were 
not calculated since the orfice area to well ratio is larger than the 


Santa Cruz well, providing response characteristics better than those 


of the Santa Cruz. A 1/4 inch orfice was used in the Santa Cruz well. 


Pale 








TABLE VI 


Computed Response Characteristics of Santa Cruz Stilling Well 
(from Robinson, 1965) 


PERTOD ORFICE DIAMETER PHASE IAG AMPLITUDE 
(in. ) (deg) REDUCTION 
(2) 
20 sec Q:25 180 oo" 
20 sec 1.00 72 45 
60 sec 0.25 75 91 
60 sec 1.00 30 5 
25 min 0.25 6 1 
25 min 1.00 0 0 
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III. ANALYSIS PROCEDURE 


A. MONTEREY 

The Monterey tidal records were digitized at 100 points per inch 
giving a sampling interval of 34.95 seconds . The discrete Fast Fourier 
Transform utilized requires 2x2" data points, and for the records 
analyzed, 2048 data points were used which gave a record length of 
71545 seconds, or approximately 20 hours. In order to remove the tidal 
influence and reduce "leakage", the records were high pass filtered by 
fitting a least-squares polynomial curve to the raw data and the curve 
was then subtracted from each ordinate point (Appendix C). An example 
of the fit to sample data points is shown in Figure 3. It was found 
that the 6th-degree polynomial provided the best fit for this particular 
tidal cycle. 

Once the tidal effects were removed, the data is subjected to the 
Fast Fourier Transform to obtain the energy density and phase estimates 
as a function of frequency. The IBM/360 standard subroutines HARM and 
RHARM were used to calculate the Fourier coefficients, from which the 
one-dimensional spectra are derived. This method and procedure is 
explained in Appendix B. 

Subroutine RHARM gives the raw Fourier'’s coefficients 


Ai 2 EB eemOny Aas Beem Geos cer ie: By = (0 


0 N 


which are then combined as 


iB. 
+ —L j = 1,2,...N 


NL > 


to give the one-sided spectral estimate 
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Figure 3 


Least Squares Curve for Santa Cruz and Monterey, 6 Dec 69 
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ii (a, - Bs) G) 


where 


T = record length (seconds) 


Z 
spectral estimate (m -sec) 


E, 
J 
The associated frequencies are 


= _ti 2 _ N-l 
FiO en eee = 


The raw spectral estimates were then smoothed over five frequencies, 
giving each new estimate 10-degrees of freedom and increased confidence 
in the spectral peaks. Smoothing over five spectral estimates removes 
spurious peaks, but does not smear the actual peaks. The 90% con- 
fidence limits for 10 degrees of freedom are interpreted as meaning 
that one can state with 90% accuracy that the actual spectral estimate 
is greater than iB, but less than SoGIDS. as the power spectrum is 


distributed according to a CHI-SQUARE distribution. 
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B. SANTA CRUZ 
The analysis procedure of the Santa Cruz data was more complex than 
the Monterey data due to two considerations: 
(1) The Santa Cruz tidal trace is recorded on curvilinear 
coordinates, but is accepted as rectilinear coordinates 
by the CAIMA 480 digitizer. 
(2) The sampling interval of Santa Cruz is 6 seconds which 
must be interpolated and matched point for point in time 
with the Monterey data, in order to compute the cross 


Spectra. 


Z2 








The unadjusted amplitudes and time increments were obtained from use 
of the CAIMA 480 digitizer, and computer program CONVERT. The analog 
recokd was recorded on Bristol Company chart 6112, shown in Figure 4, 
Padich has curvilinear coordinates (t,R9). The digitizer then recorded 
incremental pairs, in the order X and Y (of the cartesian coordinates 
of the tidal curves), every .01 inch of stylus travel. It was there- 
fore necessary to form the cartesian coordinates by summation within 
the computer. These cartesian coordinates were then converted to the 
approximately correct T, R@ coordinates by employing geometric 
relationships. 

The computation of the correct value of t. (time axis) is critical 
to the conversion. On the time scale, 0.1 inch of record was equivalent 
to a sample interval of 6-seconds. The summation of X, X. representing 


the horizontal travel was converted to the t, R@ coordinates by 


Ce oe Ro R? Bie 
rh 1 i 


az 


where, 
t. = time in curvilinear coordinates 
R = radius of curvature for arcs of constant ts 
Y. = the ordinate on the cartesian scale 


X, = time scale in equally spaced increments of t on the 
cartesian scale. 


These geometric relations are shown in Figure 5 where 


1 iz 2 


@. = TAN Y./VR - Y. 
L 1 


i 
The adjusted heights are then computed as: 
Y, = RO 
i 


where Y, is now the adjusted sea-surface elevation value. Since the 
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Space increment is .0l1 inch, the adjusted time, t. is 
1 


t. = 600 X. 
1 1 


The adjusted sea surface elevation (with reference to an arbitrary 
level) and adjusted times, were not in general equally spaced. There- 
fore, an interpolation scheme was used to obtain data points at a 
sampling interval of 34.95 seconds in oilers Ge match the Monterey data, 
point for point (Figure 6). 

A linear interpolation was per formed between time increments greater 
than 34.95 seconds. The slope of the line is 


vee Sia Sinie 


Act, - t._) 
The horizontal spatial scale for the interpolated values is, 
xX. = (i-lL)A where A = 34.95 sec. 


The resultant interpolated amplitude is, 


5 Sie 
Sige (ome isis ee 


i-l 

Effectively, the computer checks each successive time value until 
the desired interval of 34.95 seconds is exceeded, then interpolates 
between the value greater than 34.95 seconds and the prior value which 
has an increment less than 34.95 seconds. Santa Cruz data were then 
aligned with the Monterey data point for point. 

The first 2048 points of the Santa Cruz data were handled in the 
Same manner as the Monterey data. The tidal influences were removed 


by the least-squares scheme, and the adjusted data were analyzed in 


the same manner as the Monterey data. 
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Figure 6 


Linear Interpolation 


C. CROSS SPECTRA 
Given two random, stationary time-series, £,(t) and £, (t) the 


cross-sSpectra is defined as: 


co 
a =107, 
2159 (0) = a | P19 (Te dt (1) 
where, 
T = time lag 
Oo = angular frequency 
Poe, the cross-covariance between £, (t) and £, (t) given by 
T/2 
P19 (Tt) = lim { £,(t) £, (t+T) dt (2) 
T-< 
-T/2 
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where 


tT = time lag 
T =. record length 
t = time 


If the time lag tT is replaced by -T and substitution of the dummy 


variable §6 = t - T we have, 


T/2 
P19 (-T) = Lin f £,(6 +7) £, (5)d8 
-T/2 


then, the cross-covariance is an even function such that 
12 (-7T) = 51 (tT) (3) 


The cross-spectrum (1) can be written (using Euler's formula) 


e *°' = COS ot-i SIN ot 
as 
fea) [s 6) 
1 | 
3 (0) == ; 19 (T) COS otdr-i fo, SIN otdt 
=& =-&D 


2 


since P49 (7) is not, in general, an even function, therefore, the cross 
Spectrum is a complex function of o. 

The cross spectrum combines both the cosine and sine transforms of 
Py (7). From this, if the cross-spectrum can be resolved into odd and 


even complnents, it can then be written as, 


® 1 (0) = - [o even (rt) COS otdt-i f ep odd (7) SIN otdt (4) 


since the cosine (sine) transform of an odd (even) function is zero. 


For convenience we define two functions A(t) and B(t) where 


Sil 








Sigs = ag Su 
Pag aa 


Oyo (T) + Py, (7) _ 
iow! Doe 


If tT is replaced by -T in (5), and using (3) gives Ajo (-7) = “Ao (1) 
and Bio (-7) = Bio (7) where A is an odd function and B even. Since 
P19 (7) is resolved as 
Pig) Gai eta 


we see A and B are the odd and even parts respectively of 49° 


Using this equation, (4) can be written as 


1 : 
$ (0) ae | Bi (7) COS otdt-i i Ajo (7) SIN otdt. (6) 


The co-spectrum is defined as the real part of (6) 


foe) 


-i 
C159) Bie { Bio (7) SIN otTdT 


@& 


while the negative of this imaginary part defines the quadrature- 


Spectrum 
1 
Qi 9 69) mao i Aig (7) SIN pone 


The cross spectrum is then written in complex form 


If a time series is correlated with itself the resulting co- 
variance function is even. The cosine transform of the auto-covariance 
function is the power spectrum of the record. The co-spectrum is the 
cosine transform of the even part of the cross-covariance function, 


and gives the amount of power due to in-phase (0°) or out of phase + 
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180° components between the two records. The quadrature spectra is 
similar, however, the components are out of phase by a constant amount. 
The spectral phase difference E19 (0) is defined by 


ST ee 


(o) = TAN C5) 


©12 


and is the argument of the complex conjugate of the spectrum. 

In this study the Fast Fourier Transform method was used to determine 
the power spectra vice the auto-correlation method. This resulted in a 
slightly different scheme of calculating the cross-spectra, as described 
below. 


From subroutine RHARM the Fourier coefficients 


Zone 27h) eg ag ee NN 


were determined for both the Monterey and Santa Cruz data. 


From the Fourier coefficients we form the modules for each station, 


a. iL}8) : 
Fy) q Ae = IFyle “2 (Monterey) 
and 
A, iB. ; 
F.(j) = Sas |F,|e°2 (Santa Cruz), 
or 
A, iB. 
Foci) = Fo(j)e °2 = oa 5 (complex conjugate) 
where j = 0,1,2,...N. 
Then the cross-correlation function, 
19 (T) = = FoG) F,(ie FO" 
;= 


where Oo, = 2n/N dt 
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and the cross spectra 


D2 
= Re areas, 1 =1 10,7 
-T/2 


form a Fourier integral pair. 


The cross-spectra 19) can be expanded as 


5(4) = FG) FGG) = (ileal FS G)le "2 


1, (4) |, Ci) Le 1&2? 


where 
2 Tan! (-4) (Monterey) 
J 
HGS tan! (—). (Santa Cruz) 
J 


By use of Euler's formula, 
i090 e e 
|Ale = lal (cos 0+ i sin 8) 


we have, expanding 194); 


+ iF, G)1{F,G)| SIN (, - e5). 
or, eo = [co-spectra (in phase)] + i [quad-spectra (out of phase)]. 
For the co-spectra, if €, = €9: COS (O) = 1, and we have maximum 
value for in phase. If €y = eo + T1/2 in the quad spectra, SIN n/2 = 1 


and we have maximum out of phase. If we arbitrarily let Ey = O then 


es 1-00 


aa €1 ~ Cos ony CO 


and the phase difference, 


E 


eel Ly 
= TAN aor 
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From this, the co-spectra, quad-spectra, cross-spectra and phase 


difference was computed for each spectral component. 
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IV. ANALYSIS RESULTS 


A, IDENTIFICATION OF LONG WAVES PRESENT 

Spectral analysis of ten tidal records indicated that the long-wave 
activity in Monterey Bay is characterized by a few specific spectral 
peaks which dominate the spectra. Their seuied identification and 
possible earees are explained below. 

51.8 Minute Period 

According to Wilson (1965), the fundamental mode of longitudinal 
oscillation of the north-south extremities of Monterey Bay is 44.2 
minutes. If the bay is approximated as an enclosed basin of uniform 
depth and length, (Figure 7), the periods of oscillation are roughly 


calculated in the following manner 


Monterey Santa Cruz 





Figure /7 


Fundamental Longitudinal Oscillation of Monterey Bay 
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The shallow-water wave celerity is approximated by, 


C = Veh 
where, 


C = wave celerity 


acceleration due to gravity 


a> a 
t Ht 


average depth of bay 
The wave length of the fundamental period is 
L = 2w 


letting C = L/T, and substituting gives, 


2w 
a) ~ Veh 
or in general 


pr =i & = 1,2 
n nvAgh eens 


for the harmonics of the fundamental mode of oscillation. 


Substituting, 


> 


w = 20nm 


HT 


h 240 feet 


i 


g = 32.2 Ses 
gives the fundamental period of oscillation 


Ty = 46.5 minutes 


The 51.8 minute long-wave appears throughout the different spectra, 
and is evaluated as the fundamental period of longitudinal oscillation. 

36.1 Minute Period 

The fundamental mode of transverse oscillation of the east-west 
extremities of Monterey Bay is computed by Wilson [1965] to be 32.3 


minutes. Again, as in all Wilson's calculations, he assumed a nodal 
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line from Pt. Santa Cruz to Pt. Pinos. From Figure 8 the fundamental 


transverse oscillation is calculated. 


tite. eee 
aa 
a 
a 
Assumed a ar Moss Landing 
h 
Figure 8 


Fundamental Transverse Oscillation of Monterey Bay 
The fundamental wave length is 
L = 4b 


Using the shallow water wave equation 


c = Veh = = 


and substituting gives, 


For Monterey Bay, 


b = 8nm 


h = 240 ft 


the fundamental period is calculated as 








Ty, = 36.3 minutes. 


The calculated fundamental transverse mode, which is in agreement 
with! Wilson is present in Monterey four out of the five days analyzed, 
however, it only appears in Santa Cruz on 3 December 1969. 

In addition to the above, one can also consider the possible pre- 
sence of significant edge-wave activity out Snodgrass and Carrier, 
1956]. The edge-wave is a wave form in shallow water on sloping 
shelves and beaches. It is called an edge-wave because the waves 
approach with the crest perpendicular to the beach and travel along 
the coast. The waves are sinusoidal, and it is assumed the shelf 
slope is constant. The celerity of the edge-wave as a function of 


period and shelf slope is given by, 


C =f sin 8 
where, 
C = wave celerity 
T = period of ess 
8B = shelf slope 


g = acceleration due to gravity. 
Substituting C = L/T we have the fundamental period and harmonics 


of the edge-wave 





1 - fot 
n g sin gf 


(2n-1) n = 1, 2, 
where L is the wave length equal to twice the longitudinal distance 
across the bay for the fundamental mode. 


Substituting, 


g = 32.2 ft/sec 


og 








600 ft/8nm 


hole Lev) 


20nm 


the fundamental period of oscillation and associated harmonics of edge- 
waves for Monterey Bay are computed as, 
m= 35.6, 25.5, 20-c, 17-6, 15.7... menu ces 

These modes are also called "trapped eer from wave-guide theory 
as the wave energy is trapped into a continental wave guide where in 
theory the energy remains trapped indefinitely under idealized conditions. 

Although the 35.8 minute long-wave is significantly present as 
indicated above, it is not evaluated as an edge-wave as phase computations 
do not support any pronounced edge wave activity. 

27./ Minute Period 

The 27.7 minute period appears to be the most persistent long-wave 
in the bay. It appears in the spectrum having significant energy 
density on four out of five days analyzed, and is present at both Santa 
Cruz and Monterey on 6 December 1969 and 14 January 1970. The 27.7 
minute wave is also the dominant long-wave in Robinson's [1965] analyses. 
Such waves can be interpreted as siete ee [Munk, 1962]. The formula 
for the case of normal incidence on a two-step topography (Figure 9) 
consisting of the continental shelf of width A and constant slope S', 


and the continental slope of width B with slope S is 





Figure 9 


Two-Step Shelf Topography 
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_2ma 
= "3 Veh 


where, 


ar 
H 


fundamental period of the shelf wave 


constant which is dependent upon ratios S'/S and h'/h 


Le 2) 
tt 


where h'/h is the ratio of depth of break in continental 
shelf to depth of sea floor at the bottom of continental 
slope. 
g = acceleration due to gravity. 
For the case of Monterey Bay, choosing the 100 fathom curve as h' 
and average distance A equal to 8.0Onm and letting 48/n = 1.54 when 
S' << S and h' << h, the formula for the shelf-wave fundamental 


period becomes 


a, = sane = 
0 (60)1.54 vJgh' ° 

Substituting, 

A = 48000 ft (average distance to 100 fathom curve) 

g = 32.2 Fu/ieces 

h* = 600 ft 
we have, 

To = 29.6 minutes. 


Modification of A and h' in the shelf-wave formula for different 
locations around the bay changes the fundamental period only slightly. 
The remaining long-waves which are present during significant long- 
Wave activity are those of periods 22.5, 18.9, 16.3, 13.6, 10.1, and 
9.3 minutes. These are considered to be harmonics of the basic funda- 
mental period 51.8, 36.1 and 27.7 minutes, and are in general weak in 
spectral density. The significant spectral Dae found on various dates 


are summarized in Table VII. 
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B. INTERPRETATION OF INDIVIDUAL SPECTRA 

Figures 10, 12, 14, 16, and 18 show the synchronized tidal traces 
of Santa Cruz and Monterey on dates indicated. The traces are plotted 
with a sampling interval of one minute which is approximately equal to 
the Nyquist period of 70 seconds, which is the lowest period which can 
be analyzed using the 34.495 sampling rate. The tidal traces are repro- 
duced here for the first 16 hours of each period analyzed. 

6 November 1969 

The individual spectra for 6 November 1969 are plotted in Figure 11. 
The record tends to be a nei record with considerable long-wave 
activity. The weather Sond tmtone in Monterey Bay were calm winds and 
seas in the early morning, increasing to 15 kts of wind from the south- 
west about 2000. The surface air temperature was 55°F and skies were 
overcast for most of the day. The 42.5 minute non-recurring long 
wave appears to be significant (the term significant peak used in this 
study is a peak that has significant spectral density when compared 
with other peaks in the same record vice having a specific value) at 
both Monterey and Santa Cruz only on this date. Most noticable in the 
records are waves of periods 51.8, 42.5, 16.3 and 10.1 minutes which 
appear at both stations and have significant energy-density. 

3 December 1969 

Figure 13 shows the individual power spectra for Monterey and Santa 
Cruz from 0130-2130 3 December 1969. Examination of Figure 12 shows 
that significant wave activity in the bay commenced about 0950. 
Initially the surface winds were from the north at 2 kts, shifting to 
south at 15 kts by 1200. Skies were clear and surface air temperature 


increased from 45°F in the early morning to 58°F in the afternoon. The 
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Sea Surface Heights, Santa Cruz/Monterey 
0730 6 November 1969 to 0330 7 November 1969 
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Spectral Wave Analysis, 6 Nov 69 
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Sea Surface Heights, Santa Cruz/Monterey 
0130 - 2130, 3 Dec 69 
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Figure 13 


Spectral Wave Analysis, 3 Dec 69 
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spectra shows intense long-wave activity at Santa Cruz and moderate 
activity at Monterey. The effect of persistant southerly winds may 
account for this difference, however, no correlation is drawn. The 
36.1, 13.6 and 10.1 minute waves have pronounced energy at both 
locations, while the 51.8, 27.7 and 16.3 minute waves are significant 
at Santa Cruz alone. Note, that at Santa Cruz the 13.6 minute wave 
can very well be a harmonic of the 27.7 minute shelf wave while the 
16.3 and 10.1 minute waves could be a harmonics of the fundamental 
transverse mode of 36.2 minutes. 

6 December 1969 

The tidal traces for 6 December 1969, Figure 14, indicate heavy 
long-wave activity at Santa Cruz and only light to moderate wave 
action at Monterey. The winds in Monterey Bay were northwest at 16 kts 
in the early morning decreasing to 2 kts in the evening. The surface 
air temperature averaged 55°F and clear skies were observed during the 
day. Although in contrast to the situation observed on 3 December 1969 
no relation of wind to long-wave activity can be drawn during this 
period. The 27.7 minute shelf-wave is readily apparent in both records 
with the 13.6 minute harmonic present at Santa Cruz alone. The Santa 
Cruz picture is quite similar to that experienced earlier on 3 December 
1969 at Santa Cruz. Monterey is generally weak in energy density and 
the shelf-wave appears to be the only active wave of any significance. 

14 January 1970 

The analyses of wave records on 14 January 1970 are the most dis- 
tinctive of the study. The sea surface traces, Figure 16, show moderate 
to heavy wave activity at both stations. During the period analyzed, 


winds were southerly at 35 kts and heavy rain was falling in the bay. 
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Sea Surface Heights, Santa Cruz/Monterey 
2348 - 1848, 5-6 Dec 69 


49 





2 


Energy Density (meters -sec) 


2 


Energy Density (meters -sec) 


O0e 





pon | 
a 
c> 
o 
Frequency (millihertz) 
Santa Cruz, 2348 - 1848, 5-6 Dec 69 
ON 
a 





Frequency (millihertz) 
Monterey, 2348 - 1848, 5-6 Dec 69 


Figure 15 


Spectral Wave Analysis, 6 Dec 69 
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Figure 16 


Sea Surface Heights, Santa Cruz/Monterey 
0810 - 0410, 14 Jan 70 - 15 Jan 70 
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Spectral Wave Analysis, 14 Jan 70 
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Winds decreased to 2 kts in the late evening, however, surface air 
temperature remained near 50°F during the entire period. The 27.7 
minute wave is the primary long-wave present at both locations. In 
addition the 18.9 minute wave is present at Santa Cruz but extremely 
weak at Monterey. Again the relationship between weather conditions 
and long-waves present is not characterized by any definite pattern. 
The 51.8 minute wave is present at both locations, but is. weak. 

20 January 1970 

Wave activity at Santa Cruz was again more intense than that at 
Monterey. Winds were light and variable in the morning increasing to 
southerly at 15 kts with heavy rain in the early evening. The wave 
activity was similar to that experienced on 3 December 1969 showing 
moderate southerly winds, heavy wave activity at Santa Cruz and light 
activity at Monterey. However, with only two examples of this situ- 
ation, a clear pattern is not observed. The Santa Cruz spectrum shows 
some long-wave activity present at many discrete frequencies, however, 
the 51.8 and 13.6 minute periods are the only recurring wave periods 
present. In this record the 27.7 minute shelf-wave was weak at Santa 
Cruz. The 36.1 and 27.7 minute waves were observed at Monterey, how- 
ever, overall wave activity at Monterey was weak. 

Looking at the overall picture, including Robinson's [1969] data, 
it was seen that the 27.7 minute wave evaluated as a shelf-wave, was 
the dominant long-wave in the bay. This 27.7 minute period is averaged 
over 5 frequency band widths therefore appears within the range of 26.3 
to 29.4 minutes. Various values of average depth of the shelf and 
distance to the shelf, in general, yield periods which fall within this 


range. The fundamental transverse period of oscillation was observed 
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Sea Surface Heights, Santa Cruz/Monterey 
0215 - 2215, 20 Jan 70 
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Spectral Wave Analysis, 20 Jan 7/0 
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regularly at Monterey, but only once at Santa Cruz. The fundamental 
longitudinal period was significant at both locations most of the time. 
No definite relationship was drawn between atmospheric conditions and 
long-wave activity, however, a strong southerly wind seemed to lead to 
heavy wave activity in Santa Cruz. Also it was noted tnat most of the 
time, long-wave activity when present in the bay, was much more intense 


at Santa Cruz than at Monterey. 


C. INTERPRETATION OF CROSS SPECTRA AND PHASE 

Table VIII gives a summary of selected recurring long-period waves 
which have significant cross power when compared with other peaks in 
the same record. Some high energy-density long-waves which only appear 
once or twice in the entire study are excluded from the analyses. The 
cross spectra is a measure of the combined energy-density resulting 
from contributions of each power spectra but gives no information as to 
the forcing function of the wave or to whether the waves at the two 
locations are correlated with each other. Two isolated waves of the 
Same period but operating independently, each having considerable 
energy, would yield a high cross power. Phase relationships determine 
whether or not the waves appearing at two locations are coupled. 

Figures 20 through 24 show the cross spectrum and phase differences 
of each long-wave between Monterey and Santa Cruz. Both raw spectra 
were smoothed over five frequency band widths. However, considerable 
difficulties arose when attempting to smooth the phase relationships. 

The process of smoothing over band widths assumes that the data is 
the result of a stationary random process. The long-waves analyzed in 
this study may be viewed as deterministic functions and therefore what 


is considered in the record is a number of deterministic functions 
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superimposed on so called "white noise". The ideal way to smooth the 
phase information would be to smooth over ensembles, however, this could 
not be done because the long-waves are, in general, a transient phenomenon 
and usually not a stationary process. The Fast Fourier Transform (FFT) 
method of calculating the power spectrum is relatively new and, while 
the raw phase information is considered correct, the correct method to 
smooth relationships when using the FFT have not been completely worked 
out [Enochson and Otmes, 1968]. 

In this study the raw phase information was calculated utilizing 


the cross spectra as a weight function such that 


i = 1,2,3,4,5 


where 

e = smooth phase difference 

e. = raw phase 

2. = raw cross power. 
This method was adopted after investigating several alternatives. The 
calculation resulted in those raw phase estimates associated with a high 
cross power to be dominant in the smoothing process. Those relation- 
ships are shown in Figures 20 through 24 and summarized in Table IX for 
selected periods. There are reasons to doubt the accuracy of calcu- 
lating the smoothed phase using FFT in this manner which circumvents 
computing the correlation function. Hence, the smoothed phase calcu- 
lation must be viewed with reservations. 

Table VIII is a summary of significant cross spectral peaks for 


selected long-waves which recur throughout the five days analyzed. 


Figures 20 through 24 show the raw and smoothed cross spectra for the 
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DATE 





PERIOD 
(minutes ) 6 Nov 69 3 Dec 69 6 Dec 69 14 Jan 70 20 Jan 70 
51.8 x x x 
42.5 x 
36.1 x x 
| x x x 
22.5 
18.9 x x 
16.3 x 
13.6 x 
10.1 x x x 
2.3 x 
TABLE VIII 


Summary of Significant Cross Spectral Peaks in Monterey Bay 
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dates analyzed. As expected the cross spectra is dominated by the 27.7 
minute shelf-wave and its associated harmonics. The cross spectral 
summary is in general agreement with Table VII, the summary of indivi- 

| 
dual spectral peaks. The phase differences, between a particular wave 
at Santa Cruz and Monterey appear to be random, which implies the 
Monterey Canyon acts as a barrier dividing the bay into separate north- 
south basins which oscillate independently. This is in general agreement 


with Wilson [1965]. 


D. COMPARISON WITH WILSON'S THREE-DIMENSIONAL NUMERICAL MODEL 

Wilson [1965] assumed a boundary nodal line for his solution of his 
three-dimensional model of the oscillating characteristics of Monterey 
Bay. This node is a line drawn from Pt. Santa Cruz to Pt. Pinos on the 
Monterey Pensinula. This assumption tends to be the greatest objection 
to Wilson's results, although an enclosed bay will usually have a node 
across the seaward opening for particular modes of oscillation. The 
oscillating characteristics are depicted in Figures 25 through 28 where 
Figure 25(a) gives the general bathemetry of the bay. These depths 
were used as grid points in order to generate the numerical solution. 
The other figures show the increasing complex longitudinal modes of 
oscillation. The contour lines are the amplitudes of the water level 
normalized to the highest anti-node for the mode. Each complex mode 
of oscillation is drawn in a simplified version to the left of the 
figure. The nodal line assumption is particularly constraining to the 
longitudinal mode of oscillation which can be seen in Figure 25 where 
mode 1 (i.e., the fundamental) tends to oscillate transversely due to 
the nodal constraint. The lowest modes would be most affected by the 


nodal assumption and would have less constraint on higher modes. 
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The periods of oscillation of the longitudinal mode were calculated 

as, 
Ts me ae ac] a Sym /ateje 7) 25-5, 2UGee 2 Onae 

E45) el] AG mins, 
Mode 1 has a fundamental period of 44.2 minutes and indicates strong 
oscillation at Santa Cruz and weak at Monterey. Mode 2 is Similar to 
mode 1 while mode 3 gives the only Strong oscillation at Monterey and 
weak oscillations at Santa Cruz. Mode 4 is similar Co modes 1 and 2, 
but is somewhat more complex as are successive modes. 

Wilson does state that his assumed node may very well be incorrect 
and possible should be moved further Seaward. The boundary conditions 
of the model would then be changed, resulting in new oscillating 
periods and characteristics. This would result in a lengthening of 
the fundamental period. 

Wilson also computed two-dimensional modes of oscillation of 
Monterey Bay. This is the Cransverse oscillation of the east-west 
extremeties bounded at Moss Landing on the east and the assumed node on 
the west. The sequence of modal periods of oscillation were calculated 
as 

T, Ss 32a , WE4, SoS. 720), .. .ememanine 
The profiles of water-surface elevation for the transverse mode are 
Shown in Figure 29. 

Considering both the longitudinal and transverse modes of oscillation, 
Wilson notes that decreased Solution confidence is placed on each suc- 
cessive harmonic, however, the highest mode in which he has confidence 
is not stated. This increasing error in successive modal solutions can 
be decreased by Choosing a finer grid. Wilson's model probably is most 


valid for modes greater than mode 2 and less than mode 7. 
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As stated earlier, Wilson concluded that the deep Monterey Canyon 
does exert a profound influence on Monterey Bay oscillations as it 
effectively sets a barrier to free oscillations, between the north and 

| 
enth portions of the bay. The result is that wave activity acting in 
both locations would remain largely uncoupled and would give rise toa 
random phase. 

Table IX summarizes some of the results of this study with that of 
Wilson's model. It is apparent that the long-wave activity at Santa 
Cruz is more intense than that at Monterey. This tends to support the 
conclusion that the two basins act independently. Robinson [1965] 
found similar magnitudes at both stations and concluded that the 
oscillations were not uncoupled. The fundamental longitudinal mode of 
51.8 minutes was found to be weak at Monterey and stronger at Santa 
Cruz thus supporting Wilson's findings. The 20.5 minute wave corres- 
ponding to Wilson's mode 6 is in general weak at both stations. This 
is in agreement with Wilson's results. The greatest contrast between 
the results of this study and that of Wilson's is the effect and identi- 
fication of the 27.3 minute wave, evaluated as a shelf-wave here, but 
as the third mode of oscillation by Wilson. This wave, according to 
Wilson, should be strong at Monterey and weak at Santa Cruz. In general, 
this was not found to be the case. The ratios of energy-density at 
Santa Cruz to energy-density at Monterey is greater than 1 on 6 November, 
3 December and 6 December indicating the opposite of Wilson's findings. 
The last two periods studied show a ratio of .94 on 14 January 1970 and 
-49 on 20 January 1970, neither of which support Wilson's calculated 
results. It is considered that the 27.3 minute wave predicted by Wilson 


is in fact present during long-wave oscillations, but is evaluated as a 
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shelf-wave vice a harmonic of the fundamental frequency of oscillation 
of the bay. Most of the other periods predicted by Wilson were found to 
exist. However, it must be remembered that a discrete Fourier analysis 
j 
was performed and the data smoothed. This has the result of "forcing" 
some energy into discrete band widths allowing only certain periods to 
appear in the spectrum. 
As a result of this comparison the present study tends to verify 
many of Wilson's conclusions. The observations that wave amplitudes 
are generally higher at Santa Cruz and that the phase relationships tend 
to be random lead to the conclusion that oscillations of Monterey Bay 


at Santa Cruz and Monterey tend to operate independently of each other 


due to the impedence barrier formed by Monterey Canyon. 
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V. ERROR ANALYSIS 


Errors encountered during the course of the study were basically 
of three types: 

(1) Time errors can be introduced when recording raw data on an 
analog trace. 

(2) Errors may be generated in matching the Santa Cruz raw 
data to Monterey data. 

(3) The effect of the numerical methods of calculating the 
power-Spectrum and the cross spectrum can lead to mis- 
interpretation of results. 

Each of these effects are examined below in order to determine, in a 
qualitative sense, their effect on the calculated results. 

Time Induced Errors 

The tide gages at Santa Cruz and Monterey are both mechanical 
devices which must be wound at required intervals. Both were maintained 
within the required periods. The sampling interval of 34.495 seconds 
for Monterey and 6.00 seconds for Santa Cruz were computed by inches/ 
time relationships between two successive time checks on each record 
for each period analyzed. The times at Santa Cruz were recorded from 
radio station WWV (Naval Observatory Time) and are correct to the second. 
The times at Monterey were obtained from the telephone company and are 
accurate within + 5 seconds. Time checks were maintained daily for the 
entire period. The recording speed at Monterey was found to be 1.03 
+ .01 inches/hour and 6.00 + .02 inches/hour for Santa Cruz. Synchronized 
records are critical to phase computations. However, the resulting 


error in obtaining the raw data was considered insignificant. 
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Errors Generated in Aligning Records 


The Santa Cruz record was aligned with the Monterey record by 
converting the coordinate system from curvilinear coordinates to 
rectilinear coordinates. Linear interpolation was then used to match 
data points. 

The Santa Cruz data yielded a data point every six seconds giving 
high resolution in digitizing. The original 12,000 plus data points 
were reduced to 2048 points. Reduction and conversion of data in this 
manner resulted in very minor raw data errors. 

Numerical Calculation Methods 

The effect of approximating a continuous tidal trace with a step- 
function and transforming the step-function by numerical methods through 
use of a digital computer has several effects on the final analyses, two 
of which could result in misleading results of interpretation. 

When a continuous curve is approximated by a step function, a 
secondary wave form of high frequency energy is superimposed on the 
recorded wave form. To manatee this possible error the sampling rate 
was reduced and records were analyzed and checked for aliasing. It was 
found the aliased power was quite small and in no case greater than 10%. 

The second source of error is the effect of transforming the record 
into discrete values of energy-density as a function of discrete periods. 
Thus, the energy is forced into certain band widths and centered about 
the middle frequency. As a result of smoothing over five frequencies, 
the frequency interval was .00007 cps ranging in periods from 397.4 
minutes to 1.16 minutes. Comparison of raw and smooth power and cross 
spectra show that the raw spectra was not smeared due to smoothing the 
original 1025 spectral estimates as "spikes" in the raw spectra were 


still present after the smoothing process was completed. 
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VI. SUMMARY 


| 
j 
Power spectra analyses and cross spectral analyses of wave records 


in Monterey Bay indicate the presence of persistent unique periods 

with the phase differences between Monterey and Santa Cruz tending to 
be random. Coupled with this and the fact that wave energy at Santa 
Cruz is in general greater than that at Monterey, it is concluded that 
the Monterey Submarine Canyon has a profound effect on seiching motions 
within the bay. This effect appears to divide the Monterey Bay into 
north and south basins where oscillations occur essentially independent 


of each other. 


HLil 





APPENDIX A 


DIGITIZING GRAPHICAL ANALOG RECORDS USING THE CALMA COMPANY MODEL 480 

The Calma digitizer reduces analog graphical data to digital 
magnetic tape for computer processing and analysis. It consists of a 
freestanding tracing bed and a separate electronics/recorder module. 
The tracing bed is motor-driven for adjustment to the ROCCE comfortable 
digitizing position. 

To digitize analog graphical data directly on magnetic tape, the 
operator traces the analog curve with a moveable stylus/carriage assembly. 
As the curve is traced, movements of the stylus in either the X or Y 
direction causes pulsed to be generated by the direct digital pickoff 
mechanism. Incremental motions are recorded as characters on magnetic 
tape, leaving the task of summing the increments for whole-value coordi- 
nates to the final processing computer. This method of digitizing is 
both faster and more accurate than template methods. 

The digitizer reads and records data in the X and Y direction, with 
a sampling interval which can be set to 0.01, 0.02, 0.04, 0.08, or 0.15 
inches. The maximum absolute sampling error for the machine is 0.012 
inches. The output is external BCD, stored on 556 bpi, 7 channel, tape. 
The tape can be made compatible with the IBM/360 system, however, it is 
easier to use the CDC/6500 computer. 

A FORTRAN IV program wirtten for the CDC/6500 can be utilized for 
interpreting the digitized record. 

The program accomplished the following: 

(1) Converts the data from external BCD to display code. 


(2) Arranges the data in column matrices of 80 character length. 
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(3) Interrogated each character in the matrix to determine if 
the character is a 
(a) Point flag - a designator to label a specific data 
point of interest 
(b) Plus or minus travel in the X-direction. 
(c) Plus or minus travel in the Y-direction. 
(4) Sums the Y direction travel. 
(5) Records incremental travel in X direction (plus or minus). 

This allows for negative travel while digitizing. Negative travel 
is common if digitizing a curvilinear trace of if variations in data are 
closely spaced. 

(6) Prints the Y value V(N), for each designated increment in 
the X direction. 

(7) Prints the X increment value U(N), showing the incremental 
travel (plus or minus) for which the Y value was obtained. 

(8) Punch cards for V(N) and U(N) which are used as inputs for 
other ee. eee as desired. 

Headings may be entered on the tape with keyboard control, but 
CONVERT is not designed to read tape headings nor does the program have 
the capability to search for a particular set of data on the tape. Con- 
sequently, all data on the tape will be analyzed each time CONVERT is 
used. 

The following sequence of operations can be used as a guide for the 
digitizing procedure: 

(1) Mount chart on tracing bed so that horizontal axis of chart 
is aligned with X axis of digitizer. 


(2) Turn system power on. 
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(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 
(10) 


(11) 


(12) 


Load tape. 

Initiate "Load Forward" on tape module. 

Press RECORD ERROR. (If recorder is properly loaded, RECORD 
ERROR indicator will go out when push button is pressed,... 
digitizer will automatically shift to KEYBOARD MODE.) 

Enter required identifiers from keyboard if necessary. 

(Not required with CONVERT program.) 

Position stylus at beginning of curve. Lock Y-axis, and 
check alignment. 

Initiate TRACER Mode, de-activate SKIP (SKIP button out, 
light off), and unlock stylus. 

FIAG the tape once. 

Trace che: recond as desired. The first point which will not 
be computed, or printed, will be X = 0, Y = QO. The speed at 
which the stylus can be moved without destroying data 
accumulated is well in excess of the speed normally used to 
digitize the record. However, if exceeded, the RECORD ERROR 
will light and record will be deleted. 

When tracing is complete, lock the stylus, initiate FIAG 
PEDAL twice and IRG button once. 


Press rewind button on tape module and remove tape. 


The user of the Calma digitizer can reap the benefits of the instru- 
ment's speed and accuracy only if the digitizing operator understands 
and applies the system's basic rules of operation. It is strongly 
Suggested that new digitizing operators read, study and understand the 
textual instructions (CAIMA MODEL 480 DIGITIZER, Instruction and Main- 


tenance Manual) before attempting to operate the digitizer. 
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(13) The magnetic tape from the digitizing process is used as INPUT 


to program CONVERT. The approximate number of words on the 
tape should be estimated so that CONVERT may be dimensioned 


accordingly. 


In addition, if the tape contains any parity errors, or was not 


deguassed properly prior to usage, the program will only interpret a 


portion of the data or possibly not run at all. These problems have 


been encountered in the past. 


(14) The program CONVERT is included and is self-explanatory. 


(15) 


Output consists of (a) sampling rate, (b) length of record 
(sec), (c) amplitude values (inch) V(N), and (d) incremental 
travel in X-direction U(N). These may be summed in another 
program to determine proper time or spatial frame. 

It must be remembered that CONVERT is FORTRAN IV written for 
the CDC/6500 located at FNWC Monterey, and must be adapted 


for other computers. 
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APPENDIX B 


COMPUTATION OF FAST FOURIER TRANSFORM USING IBM/360 SUBROUTINE RHARM 
This subroutine is a one-dimensional Fast Fourier Transform (FFT) 
analysis based on the Cooley and Tukey [1965] algorithm. The program 


is designed to analyze I data points where 


I=2 Mo e3en4, 5. 2.209620: 

The FFT not only greatly reduces the number of calculations from 
earlier analysis schemes, but also reduces the round-off errors in the 
Fourier coefficients. The computations and round-off errors are 
reduced by log, (1)/1. 

The subroutine HARM is then called by RHARM to compute the complex 


coefficients: 


K = 0,1,2,...N-1. 


Then for K = 1,2,...55 - il, > (with the bar denoting conjugation) 


Le on l = 
rt ey pes 2 
SQ Seow oo a Sa She 
and, 
a TT 
me! ' At ig - 2x) 
C,. 5 (A, + AZ e Ze Ne 
for K = 1,2, = Tr 
- 1 _ At dt (Coe ry 
CuK 5 (A, AL e N-K 2°) 
for K = 1,2, = - 1 
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If we let Co = : (ReA, - ImA)) 
C Es (ReA. + ImA_) 
N 2 0 0 
We finally have: 
MS (2) B. = -2Im (CL), 
By = 0, Ay = 2 Re (Cy), 
AY = 2 Re (C,)> BY = 0. 


The Fourier coefficients Ay» Bo = 


0, A., Bi5-+-Ay, By = 0 are 


obtained for input Ki j = 0,1,2,...2N-1, where N = I/2, for the 


following equation. 


N-1l 
0 


1 : : 
=o S (A, COS (njK/N) + B, SIN (T1jk/N) + 


j ‘a 


2 Ay OD? 
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APPENDIX C 


LEAST SQUARES CURVE FITTING 

The scheme used to accomplish the least squares fit for a sixth- 
degree polynomial is described below. Given a set of data pairs (X. XY.) 
(i = 0,N) which can be interpreted as eee values of the coordinates 
of the points on the graph of y = f(x), assume that the unknown 
function f(x) -_ be approximated by a linear combination of suitably 
chosen functions, £) (x), f(x), sae £ (x) of the form F(x) = Ajfy (x) + 
A,£, (x) + A, f, (x) + A,f, (x) “+ A, £, (x) + A,f,(x) aa Af. (x) where the 
unknown coefficients Ag» Ay» ee Ag are independent parameters to be 
determined and the degree of least squares polynomial, M is such that 
M=6<N. The difference between the approximating function value 
F(xj) and the corresponding data value Y 5 is called the residual 


and is defined by the relation 


a He) = Xe (i = 0, N). 


The function F(x) that best approximates the given set of data in 


a least squares sense is that linear combination A E(x) + ALE) + 


0 


3 + AE. (x) of functions E(x) that produces the minimum value of 


the sum Q of the squared residuals where 
see? ee [ray - 1] 
Q= ; ss ; F(X.) Y, : 


For the case of curve-smoothing by polynomial least squares, we 
approximate the function y = f(x) over the range of data points 


(X,; Y,) ( i= 0, N). The parameters A,, A i Ae are then deter- 


O° 


mined such that 


86 





77 
er "| 
QF ee meee 
i st 
is a minimum. That is, a 6th degree polynomial curve is fitted to the 


data points ina least squares sense as defined earlier. The normal 


equations for the least-squares polynomial can be written as 


N+ 123 X. 2 x? ny ee A Tee 
et rl 0 i 
3 7 
ik Xk yee se BK A,/=| = X.Y 
i i l 
sx ?sx’ ys x2... rx’ A zs x.y, 
rt i 6 1s. a 


where; 
N = number of data pairs 
X, = time increment (seconds) 
Y; = amplitude values (inches) 


The solution of the matrix is the sixth-degree polynomial, 


7 it 2 3 4 5 6 
CRS Why to NC ae oi Se FAG GRIN RC GP ENC Ge INS 


which is then subtracted from each data point to remove tidal effects, 
acting as a high-pass filter. Choosing a sample interval of 90 


increments, 23 data pairs were utilized for the least squares fit. 
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